function theta1 = initializeTheta1(fullPHIon,nm,nn)

theta1.alpha = 0.001; 
theta1.rhor  = 0.8;
for i=1:2*nm+nn
    name = ['phi',num2str(i),num2str(i)];
    theta1.(name) = i*0.01;
end
if fullPHIon == 1
    for i=1:2*nm+nn
        for j=1:2*nm+nn
            name = ['phi',num2str(i),num2str(j)];
            if i~=j
                theta1.(name) = 0;
            end
        end
    end
end
for i=1:2*nm
    name = ['h0Q',num2str(i)];
    theta1.(name) = 0.005;
end
for i=1:2*nm+nn
   for j=1:i
        name = ['sigma',num2str(i),num2str(j)];
        theta1.(name) = 0.001;
   end
end
